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ABSTRACT 

We have modeled the spatial distribution of luminous X-ray binaries (XRBs) in Milky 
Way like spiral galaxies w ith an evolutionary population synthesis code developed by 
Hurley et al. (|2000L l2002h . In agreement with previous theoretical expectations and 
observations, we find that both high- and low-mass X-ray binaries show clear concen- 
trations towards the galactic plane and bulge. We also compare XRB distributions 
under the galactic potential with dark matter halo and the Modified Newtonian Dy- 
namics (MOND) potential, and suggest that the difference may serve as potential 
evidence to discriminate these two types of models. 

Key words: binaries: close - galaxies: individual: the Galaxy - galaxies: star-burst - 
stars: evolution - X-ray: binaries - stars: distribution 



1 INTRODUCTION 

X-ray binaries (XRBs) contain a neutron star (NS) or a black hole (BH) accreting from a normal companion star. They are 
conventionally divided into low -mass X-ray binaries (LMXBs) an d high-mass X-ray binaries (HMXBs) according to the masses 
of the optical companions (e.g. IVerbunt fc van den Heuvefl994h . In HMXBs, the evolved (super)giant companions, generally 
A/optical WMq, have strong stellar wind mass-loss to power a bright X-ray source for ~ 10 5 — 10 6 yr; while LMXBs, in 
which A'f pticai < 1.5M0, experience mass transfer via Roche-lobe overflow (RLOF) at a rate of ~ 10~ 10 — 10 _8 Mq yr -1 . 
Between them are inte rmediate-mass X-ray binaries (IMXBs) , in which companion stars' masses are in the range ~ 2 — IOMq 



I van den Heuvell 1975 ). Mass transfer in these binaries often occurs on a (sub)thermal timescale of ~ 10 4 — 10 5 yr through 



RLOF. 

Using distance estimates and angular distribution of LMXBs, Ivan Paradiis fc White! |l99ot) and Iwhite fc Van Paradiis! 

1 199rJ ) investigated the spatial distribution of NS and BH LMXBs in our Gala xy, and sugges t ed th at the compact objects 
had received a kick during the supernova (SN) explosions. More recent work bv lGrimm et al] (2002) using the RXTE data 
showed that HMXBs were concentrated towards the Galactic plane with a vertical scale height of ~ 150 pc while the vertical 
distributions of LMXBs was significantly broader with a scale height of ~ 410 pc, and the radial distribution of LMXBs 
peaked strongly at the Galactic bulge. But this sam ple suffers from some incom pleteness of the optical identifications/distance 
measurements at the large distances from the Sun (|Jonker fc Nelemansl 2004). Fortunately, today's sensitive, high-resolution 
X-ray observations allow the study of luminous XRBs in galaxies even beyond the Local Group, and make it possible to 
examine XRB populations in a wide range of galactic environments with different star formation histories. For example, 
XMM-Newton and Chandra observations of NGC 891, a nearby e dge-on spiral galaxy which is very similar in many respects 
to our own Galaxy, present a straightforward look from outside ([Temple et al.ll2005l ). The spatial distribution of luminous 
discrete point sources in this galaxy also shows clear concentrations towards the galactic plane and bulge. From the locations 
of 154 discrete non-nuclear ultraluminous X-ray sources (ULXs) identified in 82 galaxies observed with CTtanrfra. lSwartz et al. 
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(2004) found that the ULXs in their host galaxies were strongly peaked toward their galaxy centers. Statistical analysis of 
the X-ray point sources from the ROSAT HRI survey of nearby galaxies bv lLiu et al. I l|2006l ) showed that there is a significant 
concentration of ULXs towards galactic center in late-type galaxies. They also suggested that regular ULXs are likely to be a 
high- luminosity extension of the ordinary HMXB/LMXB population in late-type galaxies through luminosity function (LF) 
study. 

The spatial distribution of XRBs in a galaxy is determined by the initial kick velocity due to any asymmetry in SN 
explosion at the birth of a NS/BH, the galactic gravitational potential, and the mass transfer process in a binary. In the 
present work, we investigated the dynamical consequences of XRBs in spiral galaxies like the Milky Way in an theoretical view. 
We employed an evolutionary population synthesis (EPS) code to calculate t he expected number and luminosity distributions 
of XRBs in the galaxies. Then, following the approach of Paczvnski 1 1990l ). we calculated the spatial distribution of XRBs 
with luminosities > 10 3r ergs" 1 . For the galactic gravitational potential, we adopted both the standard cold dark matter 
(CDM) model and the Modified Newtonian Dynamics (MOND) model. The objective of this study is to present an integrated 
picture of XRB distribution in spiral galaxies under the two kinds of galactic potential models, and to explore the difference 
in the predicted spatial distributions, which could be testified by comparison with future high-resolution observations of XRB 
distribution in nearby galaxies. A rec ent related work is to use the detection of LMXBs in the Sculptor dwarf spheroidal 
galaxy to probe the dark matter halo i Dehnen fc Kind 20061 ). 

This paper is organized as follows. In §2 we describe the population synthesis method and the input physics for XRBs in 
our model. The calculated results are presented in §3. Our discussion and conclusions are in §4. 



2 MODEL 



2.1 Assumptions and input parameters 

We have used the EPS code developed by iHurlev et ail (|200Cl . I2OO2I ) to calculate the expected numbers for various types of 
XRB populations. This code incorporates evolution of single stars with binary-star interactions, such as mass transfer, mass 
accretion, common envelope (CE) evolution, SN kicks, tidal friction and angular momentum loss m echanics (i.e., mass loss, 
magnetic braking and gravitational radiation). Besides the modifications made by Liu fc O (12007!) to the original code, we 
have reduced the helium star wind strength by a factor of 0.6 according to lKiel fc Hurlevl ([20061 ) in modeling the formation 
processes of BH LMXBs. 

We assume that the host spiral galaxies are similar to our Galaxy, and adopt the cylindrical coordinate system (R, <f>, 
and z) centered at the galactic center. For stars born in the bulge, we simply assume that they were distributed uniformly 
between -R m in = and f? ma x = 2 kpc, while the star formation rate (SFR) in the disk varies exponentially with R, i.e., 
oc exp(— R/R cxp ) with -R oxp = 4.5 kpc, from 7? m m = 2 out to J? m ax = 15 kpc. 

Considering the different star formation processes in the galactic disk and bulge (described in §3.1), we have calculated 
the populations of X-ray sources in the bulge and disk separately. According to iBallero et al.l l|2007l ) we assume that stars in 
the bulge formed at the age of 0.4 Gyr with metallicity of 0.001 and an inital mass function (IMF) more skewed toward high 
mass than in the solar neighbourhood. We neglect magnetic braking effect for main-seque nce (MS) stars of mass 0.8 — 1.25M©, 



since metal-poor MS stars in this mass range do not have an outer convective zone fe.g. llvanovall2006h . For the disk we take 



a fixed SFR over t he lifetime of the ga laxy and solar metallicity. The values of other adop ted parameters are the same as 
the default ones in IHurlev et al.l (|2002l ) if not mentioned. The IMF of iKroupa et al.l (|l993l ) is taken for the primary's mass 
(Mi) distribution. For the secondary stars (M2), we assume a uniform distribution of the mass ratio M2/M1 between and 1 
and of the logarithm of the orbital separation In a. Tidal effect is considered to remove any eccentricity induced in a post-SN 
binary prior to the onset of mass-transfer. 



When a binary survives a SN explosion, it receives a velocity kick due to any asymmetry in the explosion ijLvne fc Lorimer 



199J). The kick velocity is assumed to be imparted on the newborn NS with the Maxwellian distribution 



7r<7 3 W 2a 2 ' 



(1) 



where a = 265kms _1 ( Hobbs et al.ll2005l ) or 190 kms~ 1 ( Hansen et al.l | l997). The direction of the initial velocity vector is 
chosen randomly. Together with the local circular motion in our Galaxy dBurton fc Gordonll978l ). this gives the initial velocity 
vectors ur, v$, v z . After evolving for a period, the binary will turn on X-rays and can be observed if it is luminous enough. 

In the mean time the motion of the binary can b e calculated if the gal actic potential is known. In our control model we 
adop t the Galactic gravitational p otential proposed by Johnston et al.l 1 19951 ). which consists of one Hernquist bulge ( Hernquist 
1990), one Miyamota-Nagai disk (jMivamoto fc Nagailll975l ) and one isothermal DM halo potential- For the MOND potential 



We have adopted other galactic potential suggested in IPaczvnskil jl99d ) and lBinnev fc Tremainel l|l987h . and found very small difference 
in the final results. 
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we use the Shan et al. 1 20081 ) model which applies MOND correction to a Kuzmin-Hernquist bulge-disk model to study orbits 
in axisymmetric potential. The potential can be constructed as, 

-GM 



$n(R,z) = 



yjB? + {a+\z\Y + h 



(2) 



where G is the gravitational constant, a the Kuzmin length, h the Hernquist len gth and M the total mass of the lens system. 
Here we adopt M = 1.2 x 10 U M Q , a = 4.5 kpc (jBinnev fc Tremainel Il987h ■ and h = 0.7 kpc (|Hernouistl 1 19901 ) in our 
calculations. Note that wh en h — * 0, Eq. (2) recovers to the thin Kuzmin d isk model with Newtonian potential given by 
Binnev fc Tremainel (|1987l ). and when a — 0, it becomes the Hernquist model (|Hernquistlll990h . The modified gravity is then 



9 = 9n + VffJvao (3) 

where gN is calculated from Eq. (2) and ao = 1.2 x 10 cms" ~ 2 dMilgromlll983l; lMcGaughll2004l ). For comparison we also 
performed calculations with the Kuzmin model presented bv lRead fc Moord (2005). 

Due to cylindrical symmetry of the galactic potential, two space coordi nates R and z are sufficient to describe the 
XRB distributions. We integrate the motion equations (i.e., Eqs. (19a, b) in |Paczvnskilll99Cl ) with a fourth-order Runge- 
Kutta method to calculate the trajectories of the binary systems and collect the space parameters of current XRBs. In our 
calculations, the accuracy of integral is set to be 10 -6 and controlled by the energy integral. 



2.2 X-ray luminosity and source type 

In our study XRBs are simply divided into LMXBs and H/IMXBs according to the mass of the optical companion, and 
we use the mass of the secondary, Ma, of 2Mq to separate them. Typically the donor star is a MS star but giant and 
white dwarf (WD) donors are also possible. For every accreting system, the bolometric luminosity (Lboi) is calculated 
based on the average mass accretion rate (M acc ) as Lboi = r\M k 
c is the velocity of light. For persistent XRB s, we adopt Lbol 
tor" ijRappaport. Podsiadlowski fc Pfah]||2004r ) to allow super-Eddington luminosities, and the critical Eddington luminosity 
Z/Edd — 4-KGMim p c/aT = 1.3 x 10 38 mi ergs" 1 (where <tt is the Thomson cross section, m p the proton mass, and mi the 
accretor mass in the units of solar mass). We assume rj^dd = 1 and 10 for accreting NSs and BHs, respectively. For transient 
sources, the luminosities in outbursts are taken to be a fraction (rfaut) of the Eddington luminosity. For NS systems, we assume 
rjout =0.1 and 1 if the orbital period P or b is less and longer t han 1 day, respectively; for BH systems, we adopt ?7 ou t = P or b/2 4 



acc c where t) is the efficiency for energy conversion and 
min(L bo i, ^Edd^Edd), where r?Edd is the "Begelman fac- 



hr and let the maximum peak luminosity not exceed 5Z/Edd ( Belczvnski fc Taaml 



2003 



The X-ray duty cycle (DC) is taken to be 0.01 empirically (e.g. Taam. King fc Ritterl 



Chen et alJI 19971 : iGarcia et al1l2003l ) 



2000) 



Finally a bolometric corr ection factor 7?boi is introd uced to convert the bolometric luminosity to the X-ray luminosity in 
the 2—10 keV energy range l|Belczvnski fc Taaml liooi ). Generally the correction factor is ~ 0.1 — 0.5 for different types of 
XRBs, and here we adopt ?7boi ~ 0.3 as our standard value. So we can obtain the simulated X-ray luminosity form as follows: 



Lx,2-10kev 



(4) 



??boi^outiEdd transients in outbursts 

fjbol min(I/bol, ^Edd-kEdd) persistent systems. 

To discriminate transien t (t) and persistent ( p) sou rces, we adopt the criteria of Ivan Paradiis 1 19961 ) for main sequence and 
red giant donors, and of Ivanova fc Kaloeeral 1 20061 ) for WD donors, respectively. 



3 RESULTS 



3.1 X-ray luminosity functions (XLFs) 



We first cons i der th e XRBs in our Galaxy. For disk sources, from t he SFR ~ O^SM pyr" 1 for stars massive than 5Mq in 
Grimm et all (|2003h . we derive the total SFR in the Galaxy following iLiu fc Lil l|2007l ). and get the specific SFR of the binary 
star population Sb = 1.028 yr _1 f or a binary fraction / = 0.5. For bulge sources, we construct a ph enomenological func tion 



of SFR(t) accordiing to Fig. 2 in iBallero et al] 1)20071 ). and assume the binary fraction / = 0.05 (jlvanova et al 



20051 ). In 

Fig. 1, we show the simulated cumulative XLFs of H/IMXBs (left) and LMXBs (middle) in the Galaxy with qceA = 0.15 



1 Dewi fc Ta uris 2000). Note that the X LF of H/IMXBs is significantly flatter than that of LMXBs, as indicated in the observed 



XLFs derived bv lGrimm et al.1 |2002j ). The breaks between 



W 5 



10 ergs in both XLFs are related to the maximum 



Eddington luminosities of various types of accreting sources (persistent NS H/IMXBs and transient LMXBs) calculated with 
Eq. (41 



Grimm et ail |2002h also combined the XLF of star forming galaxies in their sample - M82, Antennae, NGC 4579, NGC 



4736 and Circinus, with a completeness limit lower than 2 x 10 ergs . These galaxies have a total SFR of ~ 16MQyr , 
which exceeds the Galaxy SFR (~ O.25M yr _1 ) by a factor of ~ 65. They found that the XLFs of Galactic and SMC HMXBs 
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agree well with an extrapolation of the combined LF of the star-burst galaxies. We factitiously reset a SFR of 16.25M0yr -1 in 
our population calculations in order to examine the effect of SFR. The SFH is assumed to be the same as in the Galaxy (~ 12 
Gyr). We show the cumulative XLFs in late- type spiral galaxies of the right panel of Fig. 1. Note that H/IMXBs (dashed 
line) and LMXBs (dotted line) dominate at the relatively high (> 10 39 ergs -1 ) and low (< 10 39 ergs -1 ) luminosity in the XLF, 
respectively. 



3.2 Spatial distribution of XRBs in the CDM potential 

With the XLF constructed, we now calculate the spatial distribution of H/IMXBs and LMXBs at the age of 12 Gyr. We 
need to point out that we only consider the dynamical consequence of field binaries. XRBs formed in globular clusters have 
a different dynamical origin and are not included in this study. The results are described as follows. 

In Fig. 2, we show the radial distribution of luminous H/IMXBs (left) and LMXBs (middle) in a Milky Way like galaxy. 
As the figure shows, both H/IMXBs and LMXBs have a strong concentration in the direction to the galactic center along the 
galactic plane while there is a void of H/IMXBs in the galactic bulge. It is in rough agreement with the observed distribution 
in our Galaxy, considering that we have ignored the signatures of the Galactic spiral structure. Note that H/IMXBs are 
dominated by disk sources because sources in the bulge have all died early, while for LMXBs, both bulge (dashed line) and 
disk (dotted line) sources contribute to the population. 

In Fig. 2 we also show the radial distributions of luminous XRBs in late-type spiral galaxies with an enhanced SFR of 
16.25A/ Q yr _1 . We only include sources with luminosities > 10 38 ergs -1 . As the figure shows, the radial distributions also 
strongly cluster towards the galactic center region and declines along th e galactic plane, while there is a marked scarcity in 



the galactic bulge. These features can be compared with lLiu et al.l (|2006l ). who find that for the late- type galaxies, the radial 



distribution of detected ULXs shows a peak around 0.5i?25 (with R25 as the elliptical radius of the D25 isophote.), and the 
surface number density of ULXs decreases with radii until it flattens outside the D25 isophotes. 

In Fig. 3 is plotted the vertical distributions of H/IMXBs (dotted line) and LMXBs (solid line) in the luminosity interval 
10 37 - 10 38 ergs _1 (left panel) and 10 38 - 10 39 ergs (right panel), respect ively. Obvious l y H/I MXBs are more concentrated 
towards the galactic plane than LMXBs. This is in general accordance with Grimm et al. (2002), who suggested that the tail 



of high-2 LMXBs in the Galaxy cannot be solely due to the globular cluster component because only three out of nine sources 
at \z\ > 2 kpc are located in globular clusters. 

In the left panel of Fig. 4 we show the verti cal distributions of NS and BH LMXBs. They seem to be very similar, which 
is in general agreement with those obtained by Jonker fc Nelemansl (2004). In order to make a more accurate comparison, 



we define Rx(\z\) = Nx(zo < \z\ ^ 10)/iVx(0 < \z\ 10) to represent the ratio of the number with zo < \z\ ^ 10 kpc to 
the total number between \z\ = and \z\ — 10 kpc for a certain type (X) of XRB. In the right panel of Fig. 4, we show the 
ratio i?x(|z|) of H/IMXBs against Zq. The solid and dotted lines correspond to sources in the luminosity interval 10 37 — 10 38 
ergs -1 and 10 38 — 10 39 ergs -1 , respectively. As shown in the figure, the more luminous the XRBs are, the more concentrated 
they are towards the Galactic plane, as expected. This can be easily understood as follows. For a companion of mass M, the 
X-ray lifetime T oc M/Lx- The higher Lx, the shorter T, and the smaller displacement from the original position. 

In Fig. 5, we show the schematic side view of XRB distribution in the galaxy. The color-bar represents the normalized 
number rati o of XRBs in the R — z pla ne. Unfortunately, due to the uncertainties in determining the distance to XRBs in the 
Galaxy (e.g. Jonker fc Nelem ans 2004), it is difficult to numerica lly compare the theoretical expectations with observations of 



bright XRBs in our Galaxy. Observationally, Temple et al. ( 20051 ) presents a outside view in X-rays of the nearby edge-on spiral 



galaxy NGC 891. The spatial distribution of luminous discrete point sources in this galaxy also shows clear concentrations 
towards the galactic plane and galactic bulge. 

We have performed calculations with varied key input parameters to investigate their effect on the spatial distributions. 
The number of HMXBs does not depend on the SFH in late-type galaxies s i nce it is much longer than the total duration for 



the formation and evolution of HMXBs (< 10 7 yr). According to iLiu fc Lil (|2007f ). for young populations, binary formation 
rate and the factor of super-Eddington accretion rate allowed can affect the XLFs most prominently. Obviously H/IMXBs 
always follow the initial spatial distributions of their progenitor stars (i.e. close to the disk) due to their relatively short 
lifetimes. Another important parameter in the evolution of close binaries is the CE efficiency qce- Our calculations reveal 
that change of the value of qce does not significantly affect the number of H/IMXBs but has a strong influence on LMXB 
population. For small values of «ce the orbital motion of a low-mass companion during the spiral-in process may be unable to 
driven off the envelope of the massive BH/NS progenitor, resulting in coalescence rather a compact binary. Variations in the 
CE efficiency can change the relative numbers of various types of X-ray binaries, but the main feature of XRB distribution is 
determined by the dynamical processes during the formation of a BH/NS and the galactic potential adopted. Figure 6 shows 
the radial (left) and vertical (right) distributions of LMXBs in the Milky Way with different kick velocity dispersions. It is 
obvious that there is no strong dependence of the overall spatial distribution on the kick velocity dispersion, though it can 
affect the outcome of binary evolution and the relative numbers of various types of XRBs. 
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3.3 Spatial distribution of XRBs in MOND potential 

We made similar calculations of XRB trajectories within the MOND potential. In Fig. 7, we compare the radial distributions 
of XRBs in a Milky Way like galaxy in both MOND and CDM potentials. The solid, dotted, and thick short dashed lines 
represent those in Shan-Zhao model, Kuzmin model and the CDM model, respectively. Note that LMXBs in MOND potentials 
also show clear concentration towards the Galactic bulge as in CDM model, but they are peaked at ~ 1 — 4 kpc from the 
galactic center. This feature can be seen clearly in the schematic side view (Fig. 8) of the distributions of XRBs in the galaxy. 
In Shan-Zhao model (left panel) the distribution shows a peak ~ 2 — 3 kpc from the galactic center, while in Kuzmin model 
(right panel) it peaks at ~ 3 — 4 kpc from the galactic center. The remarkable scarcity of XRBs within the galactic bulge may 
provide interesting clues to falsity one of these models through comparing XRB distributions with predictions from the CDM 
and the MOND potentials. 

We note the difference between these models is mainly due to the potential discrepancy, especially the contribution from 
the bulge component. The CDM model contains three components, one Hernquist bulge, one Miyamota-Nagai disk and one 
isothermal DM halo potential. It is the bulge potential that leads to the peaked distribution of XRBs in the galactic center. 
The Kuzmin model has only one disk component, which results in few XRBs within the galactic bulge. In Shan-Zhao model, 
the auxilary point potential in Kuzmin model is replaced with an auxilary Hernquist potential, that's to say, this model 
contains not only Kuzmin disk potential but also Hernquist bulge component. This bulge potential causes the distribution 
peak to move inwards the bulge. We also find that in CDM model there is a more extended radial distribution for XRBs than 
in the MOND model. This is mainly caused by the potential of the DM halo, which is thought to be reside beyond 10 kpc 
from the galactic center to account for the missing mass. 



4 CONCLUDING REMARKS 



This study shows that, with current understanding of binary evolutions and galactic structure, it is possible to investigate 
both the luminosity function and spatial distribution of luminous XRBs in nearby galaxies, although the results are subject to 
many uncertainties and simplified treatments. For example, in our calculations, only primordial binaries were considered while 
in dense environments like the galactic bulge, dynamical formation channels such as tidal capture, exchang e encounters, and 
direc t collisions, may play an important role in binary formation and change the distribution of the XRBs ([Voss fc Gilfanov 
2007). Additionally, we adopted a simplified radial distribution of newborn binaries in the disk and in the bulge. The actual 
initial distribution is related to the structure of the spiral arms, the distribution and evolution of the giant HII regions 
and warm CO clouds, which are still poorly known. Finally, recent dynamical encounter of gal axies may also lead to the 



prevalence of ULXs population and make a great influence in their spatial di s tributions (e.g., Fabbiano. Zezas fc Murray 
200 ll : IWolter fc Trinchieril 120031 : iBelczvnski et all 120041 ; IColbert fc Millerl I2OO5I : iFabbiano fc White! l2006h . Although detail 



comparison between observations and theoretical predications is not available at present, rough agreement can be obtained. 
In particular, our calculations show that XRBs in CDM and MOND potentials may have distinct radial distribution around 
the galactic bulge, suggesting a new way to constrain the nature of DM and test the law of gravity. Our work motivates 
further efforts to explore the origin of the spatial distributions of luminous XRBs around the galactic center regions. 
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Figure 1. Cumulative luminosity function of H/IMXBs (left), LMXBs (middle) in the Milky Way and luminous XRBs in late-type 
spiral galaxies (right). The SFR of late-type spiral galaxies is 65 times of that of the Milky Way. Note that H/IMXBs (left) 
dominate by persistent sources, both neutron star (dashed line) and black hole (dotted line) persistent sources. In the middle panel 
the dashed and dotted line represents the disk and bulge sources, respectively. In the right panel the thick solid line is the combined 
luminosity function of both LMXBs (dotted line) and H/IMXBs (dashed line). Note that H/IMXBs dominate at the relatively high 
(> 10 3!) crgs~ x ) luminosity while LMXBs dominate at relatively low (< 10 39 crgs _1 ) luminosity in the XLF. 
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Figure 2. Radial distributions of H/IMXBs (left), LMXBs (middle) in the Milky Way. Only sources with luminosity L x > 10 37 ergs — 1 
are plotted. In the right panel shows the radial distribution of late-type spiral galaxies ( L x > 10 38 , 10 39 ergs _1 ). The origin of the 
coordinate is at the Galactic Center, and R is the distance from the Galactic center. 
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Figure 3. Vertical distributions of LMXBs (solid line) and H/IMXBs (dotted line) in the luminosity intervals log(Lx,2-iukcv) = [37, 38] 
(left) and [38, 39] (right) in the Milky Way. 
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Figure 4. Left panel: Vertical distribution of LMXBs with BH (solid line) and NS (dotted line) accretors in the Milky Way. Right panel: 
Variation of the ratio Rx(M) f° r H/IMXBs in luminosity intervals log(Lx,2-l0kev) = [37,38] (solid line) and [38,39] (dotted line) with 
zq in the Milky Way. 
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Figure 5. The spatial distribution of luminous XRBs in the Milky Way in the CDM potential. The origin of the coordinate is at the 
Galactic center, R and z are the horizontal distance from the Galactic center and vertical distance from the Galactic plane in units of 
kpc. The colorbar represents the normalized number of XRBs. 
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Figure 6. Dependence of the radial (left) and vertical (right) distribut ion on kick velocit y dispersion of LMX Bs in the Milky Way . 
The solid and dotted lines correspond to models with a = 265kms~ _1 iHobbs et alj|2005l) and a = 190kms 1 l lHansen etaU 1997), 
respectively. 



10 Z. Zuo et al. 




10 12 14 16 18 20 

R(kpc) 



Figure 7. Radial distributions of XRBs with L x > 10 37 ergs _1 in the Milky Way in MOND and CDM potentials. The solid, dotted, and 
thick short dashed lines represent the results in the Shan-Zhao model, Kuzmin model and CDM model, respectively. 





Figure 8. Same as in Figure 5, but for the MOND potential with Shan-Zhao model (left) and Kuzmin model (right). 



